Electronic structure of the muonium center as a shallow donor in ZnO 
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The electronic structure and the location of muonium centers (Mu) in single-crystalline ZnO were 
determined for the first time. Two species of Mu centers with extremely small hyperfine parameters 
have been observed below 40 K. Both Mu centers have an axial-symmetric hyperfine structure along 
with a (0001) axis, indicating that they are located at the ABo y and BC|| sites. It is inferred from 
their small ionization energy (~ 6 meV and 50 meV) and hyperfine parameters (~ 10~^ times the 
vacuum value) that these centers behave as shallow donors, strongly suggesting that hydrogen is 
one of the primary origins of n type conductivity in as-grown ZnO. 



Zinc oxide (ZnO) is one of the most promising semi- 
conductors for the next generation of electronic and opt- 
electronic devices. It has already been applied to trans- 
ducers, phosphors and varistors, due to its unique piezo- 
electric, optical, and electrical properties. In these appli- 
cations, polycrystalline material has mainly been used. 
Moreover, recent progress in single crystal growth ^ has 
opened up new possibilities, like bright blue and uv light 
emitters. Optical uv lasing has already been observed 
even at room temperature ||^ . 

For applications to optoelectrical devices, it is crucial 
to control the bulk electronic conductivity of crystalline 
ZnO. However, it is notoriously difficult to obtain intrin- 
sic ZnO, ending up with materials showing strong n type 
conductivity. In spite of more than 20 years of investi- 
gations, the origin of this unintentional carrier doping is 
still controversial. It has long been speculated that the 
dominant donor is a native defect, either oxygen vacancy, 
or zinc interstitial Unfortunately, recent theoreti- 

cal investigations have revealed that none of those native 
defects behave as shallow donors 

Recently, it was theoretically pointed out that hydro- 
gen (H), which is quite diificult to remove from the crystal 
growth environment, is an excellent candidate for such a 
shallow donor ||^. As shown in Fig. 1, ZnO crystallizes 
in the wurtzite structure corresponding to an elongated 
zinc blend structure with hexagonal symmetry around 
the [0001] axis. The lattice parameters are known by ex- 
periments as a = 0.325 nm, c/a = 1.602 and u = 0.382 
in normalized coordinates. From a first-principle calcu- 
lation, the lowest energy configurations for hydrogen are 
predicted to be at the BC^ site, with a nearly equivalent 
formation energy for the BC||, ABo^j_, and ABo,|| sites 
Experimental evidence for this scenario has been 
claimed in several reports [^-11|, where an increase in 
the conductivity was observed upon introducing H into 
ZnO. 

In this Letter we report on a determination of the elec- 
tronic structure and the location of a muonium (Mu, an 
analogue of isolated hydrogen whose proton is substi- 
tuted by a positive muon) as a shallow donor in ZnO. 



By using single-crystalline ZnO, two species of muonium 
have been clearly distinguished. The muonium center is 
readily observed in a wide variety of semiconductors af- 
ter positive muon implantation, and has been serving as 
a unique source of information on the electronic structure 
of isolated hydrogen centers [|2|. While the dynamical 
aspect (e.g., diffusion property) may be considerably dif- 
ferent between Mu and H due to the light mass of Mu 
^rup), the local electronic structure of Mu is virtually 
equivalent to that of H after a small correction due to the 
difference in the reduced mass (^ 4%). It is now well es- 
tablished in elemental and III-V compound semiconduc- 
tors that there are two stable (and metastable) sites, one 
at the center of the matrix bond (i.e., BC-site, Mu^^) 
with a large outward relaxation of the nearest-neighbor 
(nn) host atoms, and the other around the center of a 
tetrahedron cage (i.e., Td site, Mu§.). While Mu!^ has 
a large isotropic hyperfine parameter (almost the same 
order of the vacuum value, = 4463 MHz), the hy- 
perfine parameter of Mu^^ has a value about one order 
of magnitude smaller with a large unpaired spin density 
distributed on the nn host atoms. Recently, a novel muo- 
nium state having an extremely small hyperfine parame- 
ter (10"'* X Afj_) has been reported in a II- VI compound 
semiconductor, CdS jl^, suggesting that such a shallow 
Mu center (and H center as well) might be present in 
ZnO to serve as a donor. 

The experiment was performed at the Meson Science 
Laboratory (located in KEK) which provides a pulsed (50 
ns pulse width and 20 Hz repetition) beam of 100% spin- 
polarized muons with a beam energy of 4 MeV. The muon 
beam with longitudinal polarization was implanted into 
a single-crystalline wafer (40 mm diameter, 0.5mm thick- 
ness, [0001] orientation) of ZnO obtained from Eagle- 
Picher Industries, Inc. The conventional time differential 
muon spin rotation (/iSR) measurements were performed 
under a magnetic field applied in two different orienta- 
tions: one in the transverse direction (TF, B in Fig. 1) 
and the other in a tilted direction {B') with respect to the 
initial muon spin polarization P^. To obtain the tilted 
field, both transverse and longitudinal (LF) magnetic 
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fields were applied simultaneously. In the case of hyper- 
fine parameter measurements, the specimen was placed 
on a cold finger with the [0001] axis parallel with (O- 
face up). As shown in Fig. 1, the [1120] axis was set either 
perpendicular to B (Fig. la) or parallel with B (Fig. lb) 
to examine the angular dependence of the hyperfine con- 
stants. For the temperature dependence measurements 
of the muonium fraction, the [0001] axis was tilted by 45° 
to Pf, (i.e.,BZ[0001] = 45°) while B ± [1120]. 

It has been inferred from TF (=2.00 mT, 4.00 mT 
and 30.0 mT) measurements that only a single diamag- 
netic muon state is present above 40 K. The relaxation 
rate is almost independent of temperature with a rate 
of ~ 0.022(6) /is~^ for Gaussian damping, which is con- 
sistent with the dipole-dipole interaction of muons with 
^''Zn nuclei (natural abundance 4.1%). On the other 
hand, the muon spin rotation signal changes drastically 
below 40 K. A typical /iSR spectrum is shown in Fig. 2 
with fit errors, where the data were obtained under 30.0 
mT with i?Z[0001] = 45°. Fig. 3 shows the temper- 
ature dependence of the FFT spectrum, in which two 
pairs of satellite lines are seen with their position situ- 
ated symmetrically around the central line corresponding 
to the precession of diamagnetic muons (with the gyro- 
magnetic ratio 7^ = 27r X 135.53 MHz/T) at 5 K. The 
splitting of these satellites remained unchanged when the 
applied field was changed to 2.00 mT or 4.00mT. More- 
over, a nearly equivalent frequency spectrum was ob- 
served when the specimen was rotated by 90° around the 
[0001] axis (i.e., between (a) and (b) in Fig. 1). These 
results strongly suggest that two muonium centers with 
extremely small anisotropic hyperfine parameters exist 
in ZnO. The hyperfine parameters are about 10~^ times 
smaller than the vacuum value and they are symmetric 
to the [0001] axis. 

Provided that the hyperfine interaction has an axial 
symmetry, we expect two muonium precession signals for 
the high-field limits with frequencies 

z._(0)^i.o-^AK0), (1) 

u+{e)^ 1.0 + 1^1.(6), (2) 
Aiy{e) = A{e) =1 A|| cos^{9) + Aj_sm^{e) |, (3) 

where 27ri/o = ^^B with B = \B\ being the applied field 
{B > 27rA/7e, where 7e = 27r x 28.024 GHz/T is the 
gyromagnetic ratio of electron), 9 is the angle between 
B and the symmetry axis [0001], and A|| and A± are 
the hyperfine parameters parallel and normal to [0001], 
respectively. It was revealed upon preliminary analysis 
that the fitting of fiSK time spectra assuming the three 
frequency components (z/q, i^-, and v^) did not reproduce 
the time spectrum, yielding a large fraction ('^ 20 %) of 
fit errors and poor reduced (— 2.40). On the other 
hand, as suggested in Fig. 3a, fitting analysis with two 



sets of satellites including five components {i/q, i^i-, and 
i^i-)-, i = 1,2) turned out to yield a satisfactory result 
with drastically improved (— 1-52). This indicates 
that there are two species of Mu centers with respective 
fractional yields in this compound. From the spectrum 
with B ± [1120] (Fig. la), the hyperfine parameters are 
deduced to be 

Ai(90°) = I Aij_ 1= 358(4) kHz, (4) 
^2(90°) = I A2± 1= 150(4) kHz. (5) 

Combining this result with the data under a tilted field B' 
(where 9 = 54.0°, Az/i = 495(2)kHz, Av2 = 298(4)kHz), 
the rest of the hyperfine parameters are deduced as 

II I = 756(13) kHz, (6) 
^2 II I = 579(19) kHz. (7) 

As shown in TABLE I, the angular dependence of the fre- 
quencies {Avi ) calculated by the above parameters is in 
excellent agreement with the experimental observation. 

The possibility that these Mu centers have a hyperfine 
tensor with the symmetry axis parallel to BCj^ is elimi- 
nated by the fact that the observed precession frequency 
is independent of the rotation of the crystal around the 
[0001] axis by 90°. Another attempt to explain this by re- 
sorting to a sufficiently small anisotropy with the symme- 
try axis parallel to BC_l fails to account for the difference 
between A|| and A± consistently with the data. Thus, 
we conclude that there are two species of Mu centers, 
both of which have axially symmetric hyperfine struc- 
ture along with the [0001] axis. Hereafter, we denote 
these two centers as Mu/ and Mu// with the correspond- 
ing hyperfine parameters, Ai(9) and A2{9), respectively. 
The static dielectric constants in ZnO are reported to 
be 7.8(3) for perpendicular and 8.75(40) for parallel to 
the [0001] axis [jl^. The degree of obtained anisotropy 
for the muonium hyperfine tensor (^50%) is much larger 
than that of the dielectric constant (^^10%), indicating 
that the anisotropy is determined by the local electronic 
structure with the BC|| and ABq || sites (see Fig. 1) be- 
ing the most probable candidates for the sites of those 
Mu centers. Considering the magnitude of anisotropy 
in the hyperfine tensors, it would be reasonable to pre- 
sume that Mu/ is located at the ABq || site and Mu// 
at the BC|| site. Let us compare our results to a simple 
model of shallow level centers in a dielectric medium. In 
this model, the hyperfine parameter is inversely propor- 
tional to the cube of the Bohr radius (a^) of the bound 
electrons. The isotropic part of the hyperfine parameter, 
Aiso (= lA\\ + |yl_L), is 491 kHz for Mu/ and 293 kHz 
for Mu//. Compared with = 4463 MHz, one obtains 
fld = 21ao = 1.1 nm for Mu/ and ad — 25ao — 1.3 nm 
for Mu// (where oo is the Bohr radius of the free Mu). 
On the other hand, the Bohr radius for a hydrogen-like 
defect is calculated from the average dielectric constant, 
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e = 8.12, and the electron effective mass, m* = 0.318me, 
of ZnO ||l^, i.e. ad = {e/me/m*)ao = 25.5ao. This value 
is qualitatively in good accord with those of Mu/ and 

Mu/7. 

The temperature dependence of the amplitudes of 
Mu/, Mu//, and diamagnetic muon are plotted in Fig. 4. 
The total yield of all states are almost independent 
of temperature, suggesting that Mu/ and Mu// are 
ionized to a diamagnetic muon above the transition 
temperature(^ 40 K). It is unlikely that these Mu en- 
ergy levels are just above the valence band. Otherwise, 
the temperature dependence of the muonium charge state 
would not be expected due to the n type conductivity 
of the present specimen where the Fermi level is much 
higher than the mid-gap level. These results indicate 
that the Mu centers act as shallow level donors. Thus, 
since Mu centers simulate the electronic structure of H 
in ZnO, our result provides convincing evidence that the 
hydrogen centers in ZnO are shallow donors, leading to 
n type conductivity in ZnO. 

The activation energies of Mu/ and Mu// were ob- 
tained to be 3 meV and 25 meV, respectively from the 
data in Fig. 4. According to the analysis in [^, the rela- 
tion Ed — 2Ea is satisfied between the defect level energy 
{Ed) and the activation energy {Ea), which leads to the 
respective defect level energies of Mu/ and Mu// to be 
6 meV and 50 mcV. The latter is fairly consistent with 
the calculated value of the hydrogen-like impurity model, 
13.6(m*/me/e^) = 66 meV, and the observed value of 
61mcV attributed to H in an earlier report Consid- 
ering the large ambiguity in determining the defect level 
energy for another donor at 31 meV which has a much 
lower concentration in p3| , Mu/ may correspond to this 
shallower donor. 

The reason for the absence of Mu centers at other in- 
terstitial sites is yet to be understood. Another issue is 
that a large fraction of diamagnetic muons (about 50%) 
exists even at the lowest temperature. One of the pos- 
sibilities is that muon-oxygen bounding is formed, which 
has been commonly observed in various oxides. The other 
is that the diamagnetic centers may correspond to those 
at the BCj^ or AB^ sites, where their defect energy lev- 
els are in the conduction band and/or their hyperfine 
parameters are too small to observe in our experiment. 
Further experiments, including ^SR measurements at dif- 
ferent geometry, would be helpful to address these issues. 
Meanwhile, more accurate theoretical investigations are 
strongly required to unambiguously identify the observed 
Mu centers. 

In summary, we have demonstrated that two species of 
muonium centers are formed in ZnO below 40 K with ex- 
tremely small hyperfine parameters. These centers have 
an axially symmetric hyperfine interaction around the 
[0001] axis. The temperature dependence of their frac- 
tional yields indicates that they act as shallow donors, 
strongly suggesting that hydrogen is the primary origin 



of unintentional n type conductivity in ZnO. 
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Note added: After submission of this manuscript, the 
presence of a muonium state in the powder sample of 
ZnO was reported by a separate group [ p^ . 
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FIG. 1. The crystal structure of ZnO and the geometry of 
/uSR measurements. The [0001] axis is parallel with the initial 

muon polarization P^, while the transverse field B is either 
(a) perpendicular to or (b) parallel with [1120] axis. Tilted 
field B' is in a plane defined by B and Pfj,. "BC" refers to the 
bond center sites and "AB" to the anti-bonding sites. 
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FIG. 3. FFT amplitude of the time spectrum at (a) 5.0 K, 
(b) 20 K, and (c) 46 K, obtained with TF=30.0mT, where 
arrows indicate the satellite peaks. 
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FIG. 2. fiSR time spectrum in ZnO at 5.0 K, where the 
external field B {\B\ = 30.0 mT) was apphed 45° to the [0001] 
axis. A fitting result with two species of rnuoniurn centers are 
shown with corresponding errors (difference between data and 
fitted curve). 
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FIG. 4. The fractional yield of Mu/ (open triangle), Muj/ 
(open square) and a diamagnetic muon (closed circle) versus 
temperature in ZnO. 



TABLE I. Hyperfine splitting of muonium centers in ZnO, 
where the measured values are compared with those calcu- 
lated from a particular set of hyperfine parameters. 



configuration 8 Al/lcal Al^lexp Al'2cal Al'2exp 

B _L [1120] (Fig. la) 90.0° 358(4) 358(4) 150(4) 150(4) kHz 

B ± [1120] (Fig.la) 54.0° 495(7) 495(2) 298(8) 298(4) kHz 

B \\ [1120] (Fig.lb) 90.0° 358(4) 356(2) 150(4) 153(4) kHz 

BZ [1120] = 42.4° (lb) 47.6° 539(8) 541(4) 345(9) 350(4) kHz 
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